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FOREWORD 

This  report  contains  the  records  of  the  visits  to  the  U.  S.  Army- 
Engineer  Waterways  Experiment  Station  (WES)  by  Professors  H.  B.  Seed, 

15-16  September  1970,  and  R.  V.  Whitman,  26-27  October  1970,  to  discuss 
the  earthquake  resistance  of  earth  and  rock-fill  dams.  The  Office,  Chief 
of  Engineers  (OCE),  authorized  these  visits  under  its  Civil  Works  Program 
as  a part  of  Engineering  Study  5^0  entitled  "Earthquake  Resistance  of 
Earth  and  Rockfill  Dams." 

Engineers  of  the  Soils  Division,  WES,  actively  engaged  in  the  dis- 
cussions and  report  preparation  were  Messrs.  S.  J.  Johnson,  R.  W.  Cunny, 

L.  W.  Heller,  LT  J.  E.  Ahlberg,  and  SP5  W.  C.  Moss.  The  work  was  under 
the  general  supervision  of  Mr.  James  P.  Sale,  Chief,  Soils  Division.  This 
report  was  prepared  by  Mr.  Cunny  and  LT  Ahlberg. 

Director  of  WES  during  the  visits  and  the  preparation  of  this  report 
was  COL  Ernest  D.  Peixotto,  CE,  and  Technical  Director  was  Mr.  F.  R.  Brown. 
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11  February  1971 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Professor  H.  Bolton  Seed  Visit  15-16  September  1970,  Earthquake 
Discussions 


1.  Professor  H.  B.  Seed  visited  WES  on  15-16  September  to  discuss  the 
earthquake  resistance  of  earth  and  rock- fill  dams.  A list  of  those  who 
attended  the  discussions  is  given  in  Incl  1.  Professor  Seed  was  questioned 
regarding  his  views  on  earthquake  design  input,  appropriate  soil  properties 
for  earthquake  analysis,  and  earthquake  analysis  procedures.  The  remarks 
that  follow  are  the  writers’  interpretation  of  the  comments  made  by 
Professor  Seed. 


Earthquake  design  input 

2.  Geologists  are  very  important  in  selecting  the  magnitude  of  an  earth- 
quake for  a particular  site.  They  are  very  capable  in  locating  faults 
and  the  development  of  fault  history  as  well  as  presenting  the  regional 
geologic  structure.  Professor  Clarence  Allen  of  the  California  Institute 
of  Technology  and  Dr.  Lloyd  Clough  of  Woodward-Clyde  & Associates  are 
particularly  proficient  in  evaluating  the  potential  effect  of  faults 

at  a particular  site.  Geological  records,  as  opposed  to  seismological 
records,  have  an  advantage  in  determining  the  potential  activity  of  an 
area  because  the  geological  records  have  a longer  history  than  the 
latter. 

3.  For  determining  the  appropriate  rock  motion  at  a site,  Professor  Seed 
recommends  determining  the  maximum  magnitude  of  an  earthquake  on  a fault 
or  faults  likely  to  be  critical  for  the  site,  a depth  of  focus,  and 

a distance  of  the  site  from  the  fault.  The  details  of  this  procedure 
are  described  in  Professor  Seed’s  paper,  "Characteristics  of  Rock  Motions 
During  Earthquakes,"  ASCE  Journal,  Soil  Mechanics  and  Foundations  Division, 
September  1969,  Professor  Seed’s  paper,  "The  Response  of  Earth  Dams  Dicing 
Earthquakes,"  included  in  the  Proceedings  of  the  Seismic  Instrumentation 
Conference  held  in  San  Francisco  in  November  1969,  and  the  University  of 
California  Earthquake  Engineering  Research  Center  Report  on  "Rock  Motion 
Accelograms  for  High  Magnitude  Earthquakes,"  April  1969. 

4.  The  rock  motion  at  the  bottom  of  the  alluvial  valley  may  be  5 to  15 
percent  less  than  the  rock  motion  of  the  outcrops  at  the  valley  walls 
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and  depends  upon  the  amount  of  overburden  in  the  valley.  Professor  Seed 
has  used  a value  of  10  percent  for  certain  analyses  he  has  made. 

5.  Bedrock  motion  is  easier  to  predict  than  ground  motion  at  the  surface 
of  a soil  or  alluvial  deposit.  Ground  motion  can  be  very  different  for 
two  locations  that  are  near  each  other  because  of  differences  in  soil 
properties.  The  maximum  velocity  of  motion  is  limited  by  the  shear 
strength  of  the  material. 

6.  Ground  motion  amplification  can  be  determined  assuming  either  (a) 
deformable  rock  properties,  or  (b)  rigid  base  rock.  The  method  using  the 
deformable  rock  properties  was  developed  by  Kanai,  is  one-dimensional, 
and  permits  energy  to  be  radiated  into  the  rock  foundation.  The  rigid 
biise  rock  method  is  a closed-energy  system  but  can  be  used  for  either  a 
one-dimensional  or  a two-dimensional  analysis;  the  magnitude  of  the 
calculated  ground  motion  is  affected  by  the  damping  characteristics  of 
the  system. 

7.  A report  on  the  effect  of  soil  conditions  on  damage  caused  by  the 
Garacas  earthquake  can  be  found  in  a University  of  California  publication. 
Earthquake  Engineering  Research  Center  Report  No.  69-2,  entitled  "Relations 
Between  Soil  Conditions  and  Building  Damage  in  the  Caracas  Earthquake  on 
July  29,  1967." 

8.  The  following  engineers  have  worked  with  strong  motion  earthquake 
records  and  are  considered  among  the  best  qualified  for  developing  design 
ground  motion  inputs  for  specific  locations:  Housner,  Newmark,  Blume, 
Arnbraseys,  Seed,  Kanai,  Rosenblueth,  Esteva,  and  Whitman. 

9.  Three  methods  for  obtaining  an  earthquake  time  history  for  analysis 
are:  fa)  use  some  previous  record,  (b)  use  some  previous  record  which 
is  modified  by  changing  one  or  both  of  the  intensity  and  time  scales, 
depending  upon  the  existing  conditions  at  the  site,  and  (c)  generate 

an  artificial  earthquake  using  some  technique  such  as  filtered  white 
noise.  Professor  Seed  prefers  using  the  method  which  best  suits  the 
particular  site  in  question.  If  a previous  record  is  available  from 
that  exact  location,  then  that  record  may  be  adequate.  If  no  previous 
record  is  available,  then  some  type  of  record  modified  for  site  conditions 
might  be  used.  Artificial  earthquakes  which  have  been  described  by 
Jennings,  Housner,  and  Tsai  in  their  Engineering  Earthquake  Laboratory 
Report,  "Simulated  Earthquake  Motions,"  California  Institute  of  Technology, 
April  1968,  can  be  used  as  firm  ground  input  but  should  not  be  interpreted 
as  rock  motion.  Unless  carefully  constructed  by  appropriate  filtering, 
an  artificial  earthquake  motion  is  considered  least  desirable. 
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10.  The  three  most  significant  characteristics  of  earthquake  motion 
are:  (a)  maximum  acceleration,  (b)  the  predominant  period  or  periods, 
anil  (c)  duration.  Duration  is  esj)ecially  important  in  problems  involving 
soil  stability;  each  cycle  of  stress  can  cause  increases  in  pore  pressure 
which  may  cause  large  permanent  deformations. 

11.  for  major  earthquakes,  one  can  move  the  earthquake  along  the  fault 
for  its  full  length  of  break  (for  magnitude  8 earthquakes,  the  fault 
break  would  be  about  200  miles).  Ihis  movement  can  produce  very  different 
time  histories  depending  upon  where  the  fault  break  begins  with  reference 
to  the  site. 

12.  Items  which  Professor  Seed  believes  need  more  study  in  the  area  of 
earthquake  design  input  are:  (a)  effect  of  motions  with  different  pre- 
dominant periods  on  the  seismic  response  of  a structure,  (b)  the  effect 
of  different  time  histories  of  motion  having  the  same  general  character- 
istics, and  (c)  the  effect  of  the  two  horizontal  acceleration  time  histories 
which  are  similar  but  have  very  different  response  spectra. 

Soil  properties  for  earthquake  analysis 

13.  The  two  most  important  properties  necessary  for  finite  element  analysis 
are  damping  and  dynamic  shear  modulus.  Although  Poisson’s  ratio  is  also 
required,  accuracy  of  this  property  is  not  critical.  Both  damping  and 
shear  modulus  vary  as  a function  of  shear  strain  and  should  be  selected 

on  the  basis  of  estimated  shear  strain  for  analysis  purposes. 

14.  Dynamic  shear  modulus  can  best  be  estimated  from  field  in  situ  tests. 
For  some  soils  such  as  saturated  natural  clays,  disturbance  can  have  a 
great  effect  on  modulus  and  if  this  property  is  measured  in  the  labora- 
tory, it  should  be  adjusted  to  more  accurately  represent  the  in  situ 
conditions.  A field  test  used  by  Weston  Geophysical  for  determining 
shear  modulus  involves  propagation  of  a shear  wave  from  one  borehole  to 
another.  Shear  wave  velocities  can  also  be  determined  by  surface  vi- 
brators and  second  arrivals  from  refraction  seismic  tests. 

15.  Damping  cannot  be  accurately  evaluated  from  field  tests,  so  lab 
tests  must  be  run  to  determine  this  property.  Different  tests  are  needed 
to  determine  the  relationship  of  damping  over  a wide  range  of  shear  strain. 

16.  Cyclic  load  triaxial  tests  are  used  to  evaluate  the  response  of 
the  soil  to  the  repeated  earthquake  loading.  For  sands,  either  undis- 
turbed samples  or  samples  remolded  to  appropriate  densities  can  be 
used.  The  confining  pressure  and  initial  axial  stress  applied  to  the 
specimens  should  cover  a range  of  stress  conditions  dictated  by  the 
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calculated  stress  conditions  in  the  embankment,  ilie  cyclic  stress  should 
be  equivalent  to  that  produced  by  the  earthquake  loading. 

17.  Susceptibility  to  liquefaction  fur  cohesionless  iiaterial  can  be  de- 
termined with  the  cyclic  triaxial  test  by  determining  the  natural  density 
of  the  deposit  and  running  a cyclic  test  on  a specimen  prepared  at  that 
density.  A small  change  in  density  near  maximum  makes  a large  difference 
in  the  number  of  cycles  to  failure  during  the  cyclic  load  test.  Relative 
density  is  raMier  difficult  to  determine  because  different  procedures 
are  used  in  different  laboratories  and  this  results  in  different  maximum 
and  minimum  densities.  A one  Ib/cu  ft  change  in  dry  density  may  cause  a 
5 to  10  percent  difference  in  relative  density  for  some  soils. 

18.  One  difficulty  with  cyclic  triaxial  tests  is  the  90  degree  rotation 
of  the  principal  stresses  during  the  application  of  the  deviator  axial 
stress  when  initial  stress  ratio  is  either  one  or  near  one.  When  the 
confining  pressure  exceeds  the  axial  stress,  extension  and  necking  de- 
velop in  the  specimens  and  nonreprosentative  conditions  are  produced 
during  the  test,  lliis  situation  becomes  aggravated  as  pore  pressures 
are  developed  and  the  effective  confining  pressure  exceeds  the  effective 
axial  stress  by  increased  amounts  and  increasing  strains  in  the  specimen 
develop  rapidly. 

19.  Dr.  Casagrande  has  pointed  out  that  pore  pressure  concentrations 
build  up  in  the  specimen  around  the  cap  in  the  cyclic  load  triaxial  test; 
this  results  in  a conservative  value  of  strength  (failure  or  large  de- 
formations at  fewer  cycles  of  load).  Professor  Seed  also  points  out, 
however,  that  using  an  incorrect  value  for  confining  pressure  (a,  ~ °q) 

at  the  beginning  of  the  test  leads  to  an  unconservative  value  of  strength 
(an  increase  in  number  of  cycles  for  failure).  Taking  both  of  these  con- 
siderations  into  account,  the  cyclic  shear  stress  which  will  cause 
failure  on  a given  number  of  cycles  for  the  cyclic  triaxial  specimen  is 
greater  than  that  for  actual  field  conditions;  for  this  reason 
Professor  Seed  reduces  the  cyclic  shear  stress  causing  failure  by  a 
factor  of  0.6  when  the  cyclic  tests  are  run  with  a = a The  basis 

U.  o 

for  this  correction  is  described  in  a report  by  Seed  and  Idriss  entitled 
"Applicability  of  Laboratory  Test  Procedures  for  Measuring  Soil  Liquefac- 
tion Characteristics  Under  Cyclic  Loading." 

20.  The  cyclic  loading  simple  shear  test  appears  to  provide  a better 
means  for  deteimining  liquefaction  characteristics.  This  apparatus  more 
closely  simulates  field  conditions  in  that  it  utilizes  the  correct 
initial  stress  conditions  and  the  rotation  of  the  principal  plane  is 
more  like  that  in  the  field;  however,  there  are  difficulties  with 


4 


WESSD  1)  j-y  I'j/i 

SUBJECT:  Frol'essar  H.  Bolton  Seed  Visit  IS-K.  SoptoinBor  l')70,  i;artlunKik.- 
Discussion.-' 

assembling  the  apparatus  and  ruiuring  the  tests,  and  great  rare  iiuist  he 
taki'ii  in  working  with  the  .pparatns. 

Earthqn^ike  analysis  procedures 

'21.  Iliere  are  various  mi‘thuds  lor  eartluiuake  aiuilysis  availal)le  to 
determine  tlie  seismic  respt)nse  of  oartli  dams.  I’rolessor  Seed  prelerti 
to  use  the  aiiproach  of  (a)  determining  tin*  stresses  in  tlie  emlsinkjiienl 
due  to  an  earthcjuiike  loailing,  (b)  subjt'ct  ing  a lalxtratory  spei'iiiien  to 
an  eciuivalent  number  of  appropriate  cyclic  stresses,  (e)  evaluating  the 
strain  obsei-ved  from  the  laboratory  specimens,  and  r<*ial  ing  tliis  to  the 
response  of  the  embankment.  One  examj)Je  ol'  such  .m  aiuilysii;  is  g.lveii 
in  Incl  2. 

22.  .Another  method  of  analysis  is  to  construct  a circular  arc  through 
the  dam  and  compute  seismic  coel'f ic lent s at  various  times  during,  thi> 
excitation  and  from  this  cfjmpute  a minimum  factor  of  safety.  See 
Incl  3.  Professor  Seed  has  scxnet  imes  ii.'^ecl  this  method  to  male  simple 
an.'tly.'ses  of  small  stnictures. 

23.  Guidelines  or  arbitrary  criteria  that  can  he  used  for  tlie  dynamic 
aik'ilysis  are:  (a)  if  th<?  shear  wave  velocity  is  less  than  .3000  |i/sec, 
the  material  may  hi;  considered  soil;  if  it  is  more  than  .'itiOO  (t/uec,  It 
may  be  considered  rock;  (b)  a magnitudr-  7 eart  Isinali-  has  approximately 
ten  equivalent  cycles  of  loading,  Ic)  a magnituile  H eartlspiake  has  aji- 
proximately  30  eriuivalent  cycles  oi’  Jfjailing;  IiJ^  stm’ 1 in  a d;im  sub.Je<;led 
to  5 percent  strain  could  be  considt.-red  stable;  (f)  soil  in  a dam  sub- 
jected to  a 20  percent  strain  is  no  doubt  unsatisfactory;  II  ) fioll  In  a 
dam  subjected  to  7 to  10  percent  strain  in  jirobably  'iVJty,  but  thl?;  would 
depend  upon  the  particular  circumstances. 

24.  Professor  Seed  believes  that  ProfVfSsor  .V-wimirk's  method  lor  detenniir 
ing  the  deforrriation  of  slopes  frtjin  seismic  e/cilat  ion  can  be  used  with 
confidence  for  cohesionless  materials  v/here  pore  jiressures  (\<i  ntd  deve|oj». 

25.  Professor  Seed  reviewed  tlu*  proposed  earthqu/ikx/  ■itmlynii-.  f'lr  Warm 
Springs  Dam.  An  acceptable  analysis  is  shfjv/n  in  Incl  4. 

Miscellaneous  comments 

26.  Tiie  reservoir  should  have  little  effect  on  the  s‘;Isfnic  refipfau-.e  of 

an  earth  dam  because  of  the  relatively  flat  embankirKmi  fjopee. ; fcxw«-ver, 
the  reservoir  loading  i.s  important  in  dfitermiining  <*f/ectjv‘.‘ 

stress  conditions. 
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27.  Judgment  must  be  used  in  evaluating  the  drainage  conditions  in  a 
lock-filldam.  In  the  case  of  Oroville  Dam,  part  of  the  embankment  was 
considered  drained  whereas,  further  in,  the  material  was  considered 
undrained.  Sometimes  it  might  be  desirable  to  make  calculations  with 
and  without  drainage  in  areas  where  pore  pressures  could  develop  to 
determine  the  significance  of  drainage  before  drawing  a conclusion  as  to 
its  significance. 

28.  For  earthquake  resistant  design,  a freeboard  of  15  ft  has  been 
used;  this  should  not,  however,  be  used  in  lieu  of  making  a good  dynamic 
analysis.  The  freeboard  height  is  dependent  upon  potential  reservoir 
landslides,  height  of  dam,  effect  of  overtopping  on  stability  of  embank- 
ment, and  consequences  of  overtopping  on  facilities  downstream. 

29.  Other  provisions  that  should  be  made  for  earthquake  resistance  are: 
filter  zones  should  be  as  thick  as  possible,  core  should  be  as  thick  as 
possible,  riprap  should  be  used  to  protect  against  erosion,  and  embank- 
ment material  should  be  compacted  to  a high  density. 

30.  Mr.  Tom  Leps  is  doing  studies  on  the  effect  of  overtopping  on  the 
flow  of  water  through  rock  filL  Bob  Weigel  at  the  University  of  California 
is  doing  work  on  model  tests  of  sliding  reservoir  slopes  and  the  resulting 
wave  action. 

31.  Professor  Seed  agreed  to  furnish  the  following  which  have  since 
been  obtained  by  mail: 

a.  A paper  by  Clarence  Allen  given  at  an  international  AEC  conference 
in  Tokyo.  This  paper  discusses  the  importance  of  using  geologic  data  in 
determining  the  seismicity  of  an  area. 

b.  The  report  of  Professor  Seed’s  findings  from  the  Caracas  earth- 
quake study;  also,  the  Weston  Geophysical  Report  containing  the  shear 
wave  velocity  profile  of  that  area. 

c.  A copy  of  a paper  given  by  Housner  at  a Geological  Gonference  on 
Reservoir  Induced  Earthquakes  held, at  Berkeley  in  May  1969. 

d.  A card  deck  of  the  computer  program  used  to  calculate  free-field 
motion  from  bedrock  motion  input  by  a lumped  mass  analysis. 

e.  A card  deck  of  the  digitized  record  of  a magnitude  eight  earth- 
quake moving  along  a fault  and  any  report  available  in  which  this  approach 
is  discussed. 
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Professor  Seed  also  agreed  to  furnish  a working  card  deck  for  his  equiva- 
lent linear  finite  element  computer  program  which  was  obtained  during 
1,T  Ahlberg's  visit  to  Berkeley  during  the  first  week  of  December. 
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♦Part-time  attendance 


Determine  the  response  of  the  embankment  using  a computer  program 
such  as  the  equivalent  linear  finite  element  method,  and  compare 
the  computed  strains  with  those  assumed  in  step  1. 

Using  a now  estimate  for  strain,  determine  new  values  of  shear  modulus 
and  damping  and  repeat  calculations  until  the  strains  computed  are 
near  the  estimated  strains. 

5.  Test  laboratory  samples  under  the  stress  conditions  encountered  in 
the  embankment. 

6.  Evaluate  response  of  embankment  from  deformation  of  laboratory 
specimens  subjected  to  the  simulated  seismic  loadings. 


j 
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SEISMIC  ANALYSIS  II 


1.  The  embankment  response  analysis  is  about  the  same  as  that  given  in 
Seismic  Analysis  I.  The  accelerations  in  the  embankment  eire  determined. 

2.  A failure  circle  is  assumed  and  a weighted  average  seismic  coefficient 
k is  computed  and  plotted  as  a time  history  for  each  slice .or  ap- 
propriate group  of  slices.  Estimate  the  number  of  equivalent  cycles 
of  a weighted  average  k that  is  compatible  with  the  computed  k-time 
history^. 

3.  Find  the  stresses  on  the  base  of  the  slices  of  the  failure  circle 
and  subject  the  laboratory  specimens  to  these  same  stresses  at  the 
same  number  of  equivalent  cycles  determined  in  paragraph  2,  Obtain 
strain  of  specimens. 

4.  The  strains  at  various  points  along  the  embankment  are  assumed  equal 
to  specimen  strains.  If  a maximum  failure  strain  criterion  is 
assumed^  then  a factor  of  safety  of  the  slope  can  be  computed  as 

3 ratio  of  failure  criterion  strain  divided  by  average  computed 
strain. 
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Professor  Seed's  Comments  on 
Dynamic  Analysis  for  Warm  Springs  Dam 

If  a dynamic  earthquake  analysis  is  to  be  madc>  it  should  include  a com- 
plete treatment  as  follows: 

1.  Determine  design  earthquakes  by  consulting  with  Clarence  Allen  (Cal 
Tech)  or  some  other  equally  competent  engineering  seismologist. 

2.  Determine  the  shear  wave  velocity  of  the  foundation  and  embankment 
for  dynamic  analysis.  Conduct  seismic  field  tests  on  the  dam  founda- 
tion materials  and  on  the  two  test  embankments. 

3.  Determine  static  and  dynamic  stresses  induced  in  the  foundation  and 
the  embankment.  Obtain  equivalent  linear  program  from  University  of 
California  (finite  element  method).  Construct"  mesh,  assign  modulus 
and  damping  values,  apply  design  earthquakes  and  compute  stresses 
and  strains;  repeat  analysis  using  improved  modulus  and  damping 
values.  Full  reservoir  condition. 

4.  Perform  cyclic  loading  triaxial  tests  on  embankment  materials.  Sample 

density  same  as  fill.  Use  consolidation  ratios,  ^ > of  1.0  and  2.0, 

3 

three  different  confining  pressures  and  three  different  deviator 
s.tresses,  the  largest  being  adequate  to  cause  at  least  15  percent 
strain  during  the  number  of  significant  cycles  of  the  largest 
deviator  stress  induced  by  the  earthquakes-.  These  tests  would 
require  18  valid  sample  tests,  which  might  require  about  twice  this 
number  of  individual  tests. 

5.  Interpret  the  effect  of  the  computed  stress  history  on  the  embankment 
materials.  Assume  that  the  strain  induced  at  various  points  in  the 
dam  by  the  design  earthquake  is  the  same  as  the  strain  on  a cyclic 
loaded  sample  subjected  to  similar  stresses.  Assess  dam  safety  in 


terms  of  strain  and  prepare  report  of  findings, 
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MEMORANDUM  FOR  RECORD 

SUBJECT:  Discussions  with  Professor  R.  V.  Whitman  26-27  October  1970 


1.  Professor  Whitman  presented  two  lectures  on  26  October  1970.  The 
notes  of  the  first  lecture,  "Amplification,"  are  shown  in  Incl  1.  The 
"Choosing  of  a Design  Earthquake"  was  the  second  lecture  and  the  notes 
are  shown  in  Incl  2.  During  the  subsequent  discussions.  Professor  Whitman 
answered  a prepared  set  of  questions.  These  questions  along  with  briefs 
of  his  answers  are  given  in  Incl  3.  Additional  remarks  that 

Professor  Whitman  made  regarding  earthquake  studies  are  recorded  below. 

Soil  properties  for 
earthquake  analysis 

2.  The  three  methods  used  by  Whitman  to  determine  the  shear  modulus  of 
soil  in  decreasing  desirability  are: 

a.  In  situ,  crosshole  techniques  or  surface  vibratory  methods. 

b.  Laboratory. 

c.  Hardin’s  empirical  formulas. 

These  methods  are  for  modulus  values  at  low  strains.  An  initial  value^to 
be  used  in  a dynamic  analysis  at  an  assumed  strain  (1  x 10  to  5 x 10  ) 

is  determined  from  the  relationship 
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3.  Equivalent  linear  theory  does  not  rigorously "acoount  for  material 
properties.  Ramberg  and  Osgood*  have  developed  equaticTrs-l^at  simulate 
soil  behavior  quite  well.  Professor  Whitman  has  used  the  Ttrliawing 
equation  to  evaluate  a damping  value  which  also  accounts  for  th^?»4iated 
waves . 

( )g  2 1 

n = mode  number 

= total  damping  for  mode  n 
= internal  soil  damping 
Y = density  of  the  material 
c = shear  wave  velocity  of  the  material 
s = represents  overlying  layer 
r = represents  underlying  layer 

4.  Whitman  has  recognized  that  a susceptible  material  becomes  liquefied 
when  the  properties  of  that  material  plot  above  the  line  on  the  following 
type  of  graph; 


Assumption:  Horizontal  ground 
surface 


T - dynamic  shear  stress  on  a horizontal  plane 
— = vertical  effective  stress 

V 

RD  = relative  density  of  the  material 

* "A  Description  of  Stress  Strain  Curves  by  Three  Parameters,"  National 
101^^1943^°"™^^^^^  for  Aeronautics,  Technical  Note  902,  Washington,  D.  C., 
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Earthquake  analysis  procedures 

5.  Amplification  theory  is  very  sensitive  to  surface  conditions  and  may 
not  be  accurate  for  shallow  surface  deposits.  This  is  because  the  shallow 
layers  have  a period  far  from  the  fundamental  period  of  the  entire  deposit. 
As  an  example: 


The  best  mathematical  model  for  structures  A and  B would  be: 


Model  A 


40  ft  soft  soil 


Model  B 


T 


200  ft  firm  soil 


i 


6.  Any  relationships  between  magnitude  or  intensity  versus  acceleration 
are  averages  and  do  not  take  into  account  the  effect  of  duration.  The 


show  that  the  Parkfield  earthquake  had  a larger  maximum  acceleration  and 
a higher  maximum  response  spectrum  than  the  El  Centro  earthquake,  but 
the  El  Centro  was  more  damaging  duo  to  its  large  duration  of  substantial 
peaks.  Whitman  suggests  that  duration  included  with  maximum  acceleration 
and  response  spectrum  would  be  a more  true  indicator  of  earthquake  motion. 

7.  Whitman’s  procedure  for  design  of  a building  is  to  assume  a maximum 
acceleration  and  velocity  and  construct  a maximum  velocity  response 
spectrum  for  this  motion.  Artificial  earthquakes  would  be  generated  to 
coirespond  with  this  spectrum.  (Note  that  an  adjustment  of  amplitude 
and  not  frequency  is  all  that  is  necessary  to  change  the  spectrum.)  This 
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is  an  appropriate  method  for  buildings  because  its  modes  have  very  different 
frequencies. 

8.  VN^iitman's  additional  steps  for  earth  dam  design  are: 

a.  Look  at  the  input  motion  to  estimate  the  shear  strain  in  the  dam. 

b.  Find  a damping  value  from  the  relationship  of  damping  versus 
shear  strain. 


This  is  internal  damping  from  laboratory  values  and  one  has  to 
estimate  radiation  damping. 

c.  Estimate  modulus  as  in  paragraph  2. 

d.  Using  an  appropriate  time  history  for  input,  find  response  of 
structure  and  iterate  damping  and  modulus  as  a function  of  computed  shear 
strains.  An  average  strain  value  of  two-thirds  peak  strain  is  a good 
value  for  iteration. 

9.  The  accuracy  of  amplification  theory  decreases  as  the  soil  deposit 
becomes  deeper.  The  first  reason  for  this  is  that  the  input  becomes  more 
difficult  to  define.  The  second  reason  is  that  with  a deeper  soil  deposit 
the  higher  modes  become  more  important.  With  a deep  deposit,  one  should 
use  a finite  element  analysis  that  includes  radiation  damping,  such  as 
that  developed  by  John  Lysmer. 

Miscellaneous 

10.  The  following  items  are  those  that  Professor  Whitman  believes  the 
Corps  of  Engineers  should  review  more  closely: 

a.  Shaking  table  tests  done  at  the  University  of  Mexico  by  George 
Prince  on  rock-fill  dams  and  the  investigation  of  the  breaking  of  particles 
under  loading. 

b.  Experiences  in  other  countries  (i.e.,  Japan,  Portugal,  Chile, 
Mexico). 
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11.  The  following  items  were  furnished  by  Professor  Whitman  to  the  WES: 

a.  "An  Investigation  into  the  Nature  of  Microtremors  Through  Ex- 
perimental Studies  of  Seismic  Waves,"  by  Ahmed  Allam. 

b.  A computer  program  "Dynamic  Fourier  Analysis  of  Layered  Systems" 
which  uses  a one-dimensional  Fourier  transform  analysis  to  compute  the 
response  of  linear,  visco-elastic,  non-uniform  soil  deposits,  subjected 
to  a base  excitation. 

c.  MIT  Civil  Engineering  Research  Report  R70-14,  "Damping  in  Soils: 
Its  Hysteretic  Nature  and  the  Linear  Approximation"  by  R.  Doby. 

d.  MIT  Civil  Engineering  Research  Report  R69-15,  "Theoretical  Back- 
ground for  Amplification  Studies,"  by  J.  M.  Roesset  and  R.  V.  Whitman. 

e.  MIT  Civil  Engineering  Research  Report  R70-37,  "Fundamental 
Period  and  Amplification  of  Peak  Acceleration  in  Layered  Systems,"  by 
G.  A.  Madera. 

f.  MIT  Civil  Engineering  Research  Report  R68-17,  "Earthquake  Simula- 
tion Models  and  Their  Application,"  by  S.  Hou. 
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Lecture  on  Amplification 
by 

R.  V.  Whitman 


26  October  1970 


1.  The  amplification  of  earthquake  motions  from  bedrock  to  the  soil 
surface  has  been  noted.  Figure  1 shows  the  acceleration  response  spectra 
tor  two  sites  relatively  close  together  with  respect  to  their  distance 

t rom  the  epicenter.  One  site  was  underlain  by  only  stiff  soil  whereas 
the  second  site  was  underlain  by  a layer  of  soft  soil.  Ihe  peak  accelera- 
tions were  quite  different  as  well  as  the  shapes  of  their  respective  re- 
sponse spectra. 

2.  A reliable  theory  is  needed  to  explain  and  predict  amplification. 
However,  not  much  data  is  now  available  for  validation  of  amplification 
tiieory.  Complicated  building  codes  are  being  introduced  and  theory  is 
needed  for  their  substantiation.  Localized  damage  in  cities  (e.g., 

Caracas)  where  it  was  not  expected  has  promoted  a closer  look  into  ampli- 
fication theory. 

3.  A comparison  of  two  amplification  theories,  wave  propagation 
solution  (Kanai)  and  lumped  shear  beam  (Seed),  is  given  in  fig.  2.  These 
are  one-dimensional  analyses  that  consider  linear  viscoelastic  material 
which,  although  it  does  not  depict  actual  soil  behavior,  simulates  the 
behavior  satisfactorily.  A necessary  assumption  for  both  theories  is 
that  horizontal  wave  fronts  propagate  vertically  to  the  free  surface. 

4.  'fhe  lumped  shear  beam  analysis,  with  the  aid  of  a high  speed 
digital  computer,  can  be  carried  out  by  using  mode  superposition  tech- 
niques or  the  more  time  consuming  step-by-step  procedures.  The  soil 
properties  of  shear  modulus  and  damping,  as  a function  of  shear  strain, 
are  necessary  for  the  computations.  A rigid  base  is  assumed  at  the  rock 
surface  which  does  not  account  for  radiation  damping.  The  mode  super- 
position technique  uses  one  value  of  damping  for  the  system.  This  means 
that  an  average  value  must  be  determined  from  the  various  soil  layers 

and  each  mode.  Most  of  the  computer  time  is  used  to  compute  the  eigenvalues, 
and  time  histories  at  any  level  can  be  produced  with  little  additional 
effort. 


5.  The  wave  propagation  theory  allows  energy  to  be  radiated  from 
the  soil  profile  through  the  bedrock.  The  effect  of  using  radiation 
is  shown  in  fig.  3.  The  dashed  curve  represents  the  velocity'  spectninr 
obtained  from  a lumped  shear  beam  analysis.  The  solid  line  represents 
the  values  obtained  from  the  wave  propagation  analysis.  The  difference 
has  been  interpreted  as  due  to  energy  being  trapped  in  the  system.  The 
soil  in  this  example  was  100  ft  thick  overlying  bedrock  and  both  methods 
used  the  same  damping  value.  In  the  wave  propagation  analysis,  a Fourier 
spectrum  has  to  be  generated  for  each  soil  layer  and  this  method  is  not 
suited  to  represent  a large  number  of  soil  layers. 


FfOCMOMi  FifiS  MJUK 


6.  'rtie  present  shortcomings  of  the  wave  propagation  amplification 
theory  are  shown  in  fig.  4.  Of  special  note  is  the  lack  of  confidence 
in  selecting  the  input  earthquake  to  be  used  with  the  theory. 

7.  Presently,  work  is  being  done  at  MIT  with  the  wave  propagation 
theory.  An  amplification  spectrum  cai  be  produced  from  the  Fourier 
spectra  produced  from  the  rock  and  soil  acceleration  time  histories 
(fig.  5).  Ihis  amplification  response  spectrum,  fig.  6,  shows  the 
natural  frequency  of  the  deposit  and  the  amount  of  damping  in  the  soil 
layer.  Ihis  damping  can  be  computed  using  three  techniques; 

a.  Amplitude  of  peaks. 

b.  Band  width  of  spikes. 

c.  Q-theory  (area  under  amplification  curves). 

Figure  7 shows  a comparison  of  damping  computed  from  data  in  Mexico  City. 
Note  that  damping  calculated  from  the  Q-theory  seems  to  agree  best  with 
the  laboratory  value. 

8.  Whitman's  concluding  remarks  are  given  in  fig.  8. 
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Lecture  on  the 

Choosing  of  a Design  Earthquake 
by 

R.  V.  Whitman 
26  October  1970 


1.  The  choosing  of  a design  earthquake  is  a complex  problem  and  it 
should  be  done  in  a combined  effort  by  a panel.  This  panel  should  be 
comprised  of  members  from  the  following  disciplines: 

a.  Seismology. 

b.  Geology. 

c.  Structural  Engineering. 

d.  Soil  Mechanics. 

2.  One  needs  to  use  a rational  approach  in  choosing  the  design 
earthquake.  Throughout  the  world  there  are  differences  in  seismicity 
and  a single  design  earthquake  cannot  be  used.  The  overdesign  of  nuclear 
power  plants  leads  to  substantial  monetary  penalties.  The  increased  con- 
struction cost  for  an  earthquake  increased  from  0.1  to  0.2  g accelera- 
tion is  one-half  to  one  and  one-half  million  dollars.  The  engineering 
design  costs  alone  are  one-quarter  million  dollars  for  an  earthquake 
analys is . 

3.  An  example  of  the  fast  rate  that  thinking  has  changed  concerning 
maximum  design  earthquakes  is  shown  in  fig.  1.  The  Parkfield  earthquake 
was  larger  than  the  maximum  probable  estimated  only  two  years  earlier. 
Since  then  even  larger  earthquakes  have  occurred. 

4.  Figure  2 shows  three  questions  which  arise  while  choosing  a 
design  earthquake.  Two  levels  of  risk,  the  operational  basis  earthquake 
(QBE)  and  the  design  basis  earthquake  (DBE),  are  presently  being  used 
and  are  explained  in  fig.  3.  Some  designers  require  a time  history 
while  others  need  a response  spectrum  (fig.  4).  Difficulties  arise  in 
using  either  input.  The  use  of  the  response  spectra  restricts  the  analy- 
sis to  mode  superposition  techniques.  A time  history  input  may  not 
include  adequate  representation  of  frequencies  most  critical  for  struc- 
tural response.  Some  firms  use  more  than  one  time  history  which  when 
combined  gives  a smoother  response  spectrum.  This  eliminates  the  peaks 
and  valleys  of  the  response  spectra  curves.  The  advantage  of  using 
artificial  time  histories  versus  an  actual  time  history  is  that  they  have 
a smoother  response  spectrum.  The  third  question  which  arises  is  where 
should  the  earthquake  be  placed  with  regard  to  the  profile.  In  the  pro- 
file of  fig.  5,  three  possible  locations  exist  for  input.  Location  1, 

in  the  bedrock,  would  give  the  best  simulation  of  soil  behavior  as  well 
as  soil-structure  interaction. 
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5.  Miignitude  is  a measure  of  tlie  size  of  an  earthquake.  Hie  more 
common  Richter  magnitude  is  determined  by  estimation  of  the  motion  of 

a standard  seismometer,  100  kilometers  Irom  the  epicenter.  Dviring  a large 
earthquake,  energy  is  released  along  the  fault  break.  Because  of  this  the 
distance  from  the  site  to  the  fault  is  more  important  than  the  distance 
Irom  the  site  to  the  epicenter  (fig.  6).  Intensity  is  the  qualitative 
measurement  of  an  earthquake  at  a particular  location.  The  modified 
Mercalli  intensity  was  originated  before  strong  motion  instruments  were 
developed  and  is  based  on  people’s  reactions  and  damage  caused  by  the 
seismic  disturbance.  Figure  7 shows  a relationship  between  intensity 
and  peak  accelerations.  The  short  dashed  line  was  proposed  by  Guttenberg 
and  Richter  and  the  long  dashed  line  is  the  more  recent  prediction  of 
Bcrshberger.  The  solid  vertical  lines  represent  data  of  some  30  earth- 
quakes for  which  measurements  of  both  intensity  and  peak  accelerations 
were  known.  This  shows  that  tiiere  is  no  good  relationship  between  the 
quantities  presented  and  care  should  be  taken  when  trying  to  predict 
quantitative  maximum  accelerations  from  qualitative  intensities.  A more 
useful  intensity  description  would  include  three  additional  criteria: 

a.  Maximum  acceleration. 

b.  Duration. 

c.  Nature  of  building  damaged. 

6.  The  AEG  presently  follows  a general  procedure  to  produce  a 
des ign  earthquake : 

a.  Find  intensity  from  historical  records. 

b.  Relate  maximum  acceleration  to  intensity. 

c.  Find  time  history  or  response  spectrum. 

In  the  areas  where  active  faults  are  present  (Galifornia),  faults  near 
the  site  are  located  and  a maximum  historical  earthquake  is  moved  along 
the  fault  to  the  point  nearest  the  site.  Empirical  charts  are  then  used 
to  determine  the  decrease  of  intensity  with  distance  to  the  site.  In 
less  active  areas  (Eastern  United  States),  a seismo-tectonic  approach  is 
used.  Ihe  maximum  intensity  is  determined  for  the  region  of  like  geology 
or  tectonics.  This  earthquake  intensity  is  then  considered  to  occur 
under  the  site  and  is  additionally  increased  one  unit  to  take  into 
account  the  possibility  of  an  even  larger  earthquake  occurring  at  the 
site.  From  this  a maximum  acceleration  is  determined.  A response 
spectrum,  such  as  Nuwmark  s,  is  then  developed  showing  the  relation 
of  maximum  acceleration,  velocity,  and  displacement  to  frequency. 

It  should  be  noted  that,  when  using  Newmark’s  chart,  high  frequencies 
(2-5  cps)  correspond  to  maximum  acceleration,  middle  frequencies 
(1/4  — 2 cps)  correspond  to  maximum  velocity,  and  low  frequencies  (less 
than  1/4  cps)  correspond  to  maximum  displacement. 
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7.  A new  approach  suggested  by  Whitman  is  based  upon  the  occurrence 
of  earthquakes  in  areas  for  which  a long  period  of  record  is  available 

and  a plot  of  intensity  versus  return  period  is  found  to  have  a constant 
slope.  By  using  an  equation  as  in  fig.  8 and  assuming  a probability  of 
failure  and  the  probability  that  the  earthquake  will  cause  failure  of 
parts  and  machines  one  can  calculate  a return  period  which  would  then 
give  a design  intensity. 

8.  A design  example  given  by  Professor  Whitman  involved  a nuclear 
power  plant  to  be  built  in  a valley  of  deep  sediments  4 to  5 kilometers 
thick  (fig.  9).  The  known  faults  are  shown  in  fig.  10  and  can  be  grouped 
into  three  systems  as  in  fig.  11.  The  problem  was  to  find  the  design 
earthquake  and  three  independent  techniques  were  used.  Figure  12  gives 
the  values  obtained  from  using  the  seismo-tectonic  approach  for  the  earth- 
quakes in  the  deep  valley  sediments.  The  maximum  historical  intensity 
was  increased  by  one  magnitude  and  this  maximum  credible  intensity  gave 

a corresponding  0.2  g maximum  acceleration.  A second  approach  was  to 
consider  the  effect  of  nearby  earthquakes  and  distant  earthquakes  (fig.  13). 
In  the  first  case,  it  was  assumed  that  earthquakes  beneath  the  deep  valley 
sediments  which  have  occurred  near  the  site  could  occur  at  the  site. 

The  effects  of  nearby  earthquakes  is  shown  in  fig.  14  and  based  upon  these 
data  and  judgment  it  was  concluded  that  a magnitude  6 earthquake  at  the 
site  would  produce  a maximum  acceleration  0.22  g and  a maximum  velocity 
7 in. /sec.  These  are  less  than  the  Parkfield  values  but  no  rupture  has 
occurred  at  the  site  in  question  and  no  intensity  that  large  is  expected. 

The  effect  of  the  magnitude  7 earthquake  located  40  km  away  was  then  con- 
sidered which  resulted  in  a maximum  acceleration  of  0.10  g and  a maximum 
velocity  of  11.0  in. /sec  as  in  fig.  IS.  Esteva’s  equations  were  used  to 
predict  the  effect  of  the  earthquake  at  40  km  distance  and  these  were 
modified  for  local  site  conditions.  Return  periods  were  used  as  the 
third  approach;  2000  years  of  historical  records  were  available.  Data 
were  plotted  for  the  last  100  years  of  data  and  also  for  2000  years  of 
record  (fig.  16).  The  thought  was  that  the  record  for  the  last  100 
years  was  the  most  accurate  and  an  intensity  of  7 or  8 (10,000  year 
return  period)  was  chosen  for  design.  The  three  approaches  are  summarized 
in  fig.  17  and  give  consistent  results.  Figure  18  shows  the  recommended 
envelope  for  response  spectrum  obtained  from  the  combination  of  the  nearby 
and  distant  fault  system.  The  maximum  acceleration  of  .2  g and  maximum 
velocity  of  20  in. /sec  were  used  as  design  values. 

9.  The  concluding  remarks  are  shown  in  fig.  19.  In  this  country, 
due  to  our  abundant  resources,  we  find  ourselves  overdesigning.  Other 
countries,  with  limited  resources,  are  designing  on  a more  rational 
basis  and  perhaps  this  line  of  thought  should  be  taken  up  in  the  United 
States. 


Incl  2 
Sheet  3 


^ ■^' 


I 


I 


c 

0) 

E e 

n Q) 

E u C 

*r-  ITS  •#- 


X I— 
ITS  CL 
Sr  VI 

♦r- 

T3 


CVJ 


CVJ 


IT) 


3 -I-  <U 

E U (/> 
•r-  O ***~^ 


X I—  C 
03  OJ  T- 
."E  > 


O 

CO 

o 

+-> 

C>si 


o 

CO 


c 

o 


E +J 
3 03 


^ ^ J 

•r-  <U  CJ( 


X 

03 


O) 

U 

o 

03 


CVJ 

o 

CO 

m 

lO 

♦ 

• 

• 

o 

o 

o 

03 

03 

•r“ 

•r* 

C 

c 

s. 

o 

o 

«f- 

tf- 

•r“ 

•r» 

1—  #0-^ 

03 

03 

03  in 

O 

0 CO 

03 

0^ 

3 

c 

C r- 

cr 

•1“ 

•r“ 

£ 

C 

•M 

T5  m 

03  -r- 

1. 

OJ  CO 

03 

•o  o> 

^ 03 

03 

S.  r— 

03  4-» 

o 

^ 03 

"O 

O 4- 

9 £ 

03  O 

S—  ‘r" 

03 

s» 

aL4-> 

•#-  CO 

V) 

V) 

M- 

• 03 

• 03 

^ r- 

X>— ' 

03 

03 

03 

s: 

z: 

O. 

32 


Incl  2 
Sheet  U 


Fig.  1 


00 


•♦-> 

• M 

3 

4J 

O 

C 

OJ  >» 

c 

■a  r- 

? 

•r-  OJ 

o 

O M- 

t- 

o ns 

<d  CO 

4-> 

CO 

OJ  c 

4-> 

<U 

01  ^ 

c 

3 o 

ns 

ns  "O 

4J 

U 

a. 

O 

c 

4-> 

C 

o 

4->  3 

s. 

•r* 

O JZ 

ns 

*o 

4-> 

c <o 

(U 

r— 

ns 

3 

S- 

4->  4-> 

o 

O 

OS 

CO  CO 

3 

Q. 

3 3 

oo 

O 

s:  E 

<u 

o 

c 


CO 

CJ 

cn 

•r— 

CO 

OJ 

c. 

s- 

CJ1 

CO 

OJ 

c 

•r>- 

o. 

*3 

D- 

•r“ 

ns 

ns 

X 

L. 

Q. 

f— • 

•r* 

<u 

•r* 

O 

3 

3 

S. 

O 

CO 

4J 

O' 

O 

o 

<U 

OJ 

c 

> 

4-> 

•r» 

O 

ns 

T? 

<v 

CO 

c 

c 

oi 

C 

•r“ 

3 

ns 

OJ 

4-) 

■o 

O 

E 

o. 

CO 

x: 

ns 

X 

3 

o 

00  “O 

OJ 

s 

(/) 


x:  (/) 

4^  (I) 

O)  I— 

(TS 
<U  C 

•—  o 

^ •r* 
fO  4-> 

^ rO 
O i- 

Sm  01  OJ 
CL  Q.^ 
O (9 
B 3 

2 o* 

E 01  ^ 

•r-  -id  4-> 
X ns  V. 
^ 3 
S O"  Oi 


nJ  w 

<U  T- 
*A 
o fd 
-Q 
-O 

•r-  C 
TJ  O) 

OJ  •r- 

</)  0> 
U 0)  ^ 
•o  rd 
E 3 

3 .•  Cy 

£ <u  x: 

•r-  ^ 4-> 

X Id  s- 

ns  3 Id 
s:  CT  o; 


31+ 


Fig.  3 


Incl  2 
Sheet  6 


J iU 
< 0 
0 
^ iU 
I'i 


j 


10  /C>/ 

r/\7, 

'C^// 

/^// 

r 


Incl  2 
Sheet  14 


;?  0/ 

K g 10 

I of  >n 
5 S o 

3 3 1/) 

s 2, 

J Ul 

5t  p 

X 2 
(P  ^ 


>1  ^^ 

s s 

O'  — 


v-*  7 /> : 

**  * "7  /x''-'  ' 

X x\7  /\'x 

XX 

A^  X'  X ' 

/ 4X'X\ 

^ / Xx"'  o"’ 

^ jO/  X' s ^ 

xx-C  / r:*  -v  « 

/ /^ : u.  \ 
/ 4'^  o-  ^ s 

.C  K u K ' / 

A 0 n '^V 

k /\  I A \/ 

!v'  •J' 


Fig.  11 


I 


o 


o 


3X15  JLV 

1 Qa»rJ)aov4 


t/? 

& 

tu 

> 

1 

0 

V 

lU 

d 

2 
fJ 

D 

h 

Ui 


Fig.  16 


>*T 


Incl  2 
Sheet  19 


Questions  and  Answers 
From 

R.  V.  Whitman  Discussions 
26-27  October  1970 


Desitiii  Earthquake 


1.  I'or  analysis,  is  an  artificial  or  an  actual  earthquake  the  best? 

One  artificial  record  can  simulate  many  actual  records.  Allen  Cornell, 
MIT,  is  capable  of  easily  producing  the  artificial  records. 

2.  How  many  earthquake  inputs  should  be  used? 

Three  or  j'uir,  as  a minimum. 

3.  If  more  than  one  input  earthquake  is  to  be  used,  should  these  be  varied 
according  to  duration  and/or  amplitude? 

Amplitude,  if  one  is  trying  to  produce  a smooth  response  spectrum. 

4 . How  does  one  .judge  when  enougli  eartliquakes  have  been  used  to  analyze  a 
structure? 

When  one  has  a smooth  response  spectrum. 

5.  What  is  the  validity  of  scaling  the  L’l  Centro  earthquake  for  designs 
in  the  Midwest?  Isn't  this  a common  practice? 

This  is  common  practice  although  it  has  no  validity. 

6.  How  does  one  develop  an  adequate  time  history  from  a response  spectrum? 
Allen  Cornell,  MIT,  is  very  proficient  at  this. 

7.  With  respect  to  the  many  artificial  earthquakes  proposed,  has  anyone 
made  an  engineering  analysis  of  the  input  evidence  (data  and  assump- 
tions) to  evaluate  the  validity  of  various  artificial  earthquakes? 

Not  to  Whitman's  knowledge. 

8.  What  is  your  opinion  of  using  ID  lumped-mass  analysis  for  obtaining 
soil  layer  response? 

This  is  a good  tool  for  iterative  purposes  but  does  not  take  into 
account  radiation  damping. 
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' Design  Analysis 


10.  What  is  the  validity  of  using  Ambraseys’  k values  from  the  envelope 
of  numerous  earthquakes? 

This  approach  does  not  take  into  account  particular  site  conditions 
and  foundation  effects.  A shear  wedge  analysis  may  be  sufficient 
when  using  this  approach  in  lieu  of  the  finite  element  method. 

11.  What  is  the  validity  of  using  k factors  in  a dynamic  analysis? 

This  is  not  a desirable  approach. 

12.  How  does  Newraark’s  method  for  determining  deformations  compare  with 
methods  proposed  by  other  people? 

Whitman  has  used  this  method  but  has  made  no  particular  comparisons. 

13.  What  are  the  "keys"  to  the  most  critical  times  during  an  earthquake? 


a. 

/o  plastic  ? 

d. 

d max? 

b. 

g max? 

e. 

stress  amplitude? 

c. 

V max  ? 

f. 

strain  amplitude? 

No 

comnent. 

14.  What  determines  a reasonable  "cutoff"  time  for  a dynamic  analysis? 
When  no  further  change  occurs  in  the  response  spectrum. 

15.  What  amount  of  permanent  deformation  or  strain  is  excessive  for  a 
plane  strain  finite  element  program? 

Not  known. 

16.  For  dams,  what  is  the  most  important  condition  for  analysis  (e.g., 
after  construction,  steady  seepage,  or  rapid  drawdown)?  Can  one 
tell  before  running  the  analysis? 

This  depends  upon  the  particular  investigation. 

17.  What  is  the  effect  of  the  reservoir  on  the  stability  of  tbe  dam? 
Do  you  know  of  any  work  being  done  in  this  area? 

Not  known. 

No. 
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13. 


What  interpretational  procedures  can  be  used  to  relate  seismic  re- 
sponse observed  from  small  earthquakes  to  anticipate  response  under 
larger  earthquakes? 

No  comment. 

19.  What  are  some  recommendations  regarding  defensive  design? 

a . Freeboard  requirement? 

b.  Thickness  of  filter  zones? 

c.  Thickness  of  core? 

d.  Special  provisions  for  spillway  and  outlet  works? 

e.  Riprap? 

The  usual  earthquake  design  provisions  should  be  made  and  also  non- 
erodible  materials  in  zones  where  cracking  is  anticipated  and  on  the 
downstream  ftice  should  be  used. 

20.  What  are  your  comments  on  evaluation  of  landslide  stability  and 
feasible  methods  for  estimating  effects  of  potential  wave  action? 

No  comment. 

21.  What,  do  you  feel,  are  the  most  critical  structures  for  a dam? 

No  comment. 

22.  How  much  design  effort  is  reasonable  for  the  earthquake  problem 
as  compared  to  a static  analysis? 

The  first  few  times  a dynamic  analysis  is  made  it  may  be  very 
costly,  but  with  experience,  this  cost  should  decrease. 

23.  What  information  do  you  have  on  reservoir  induced  earthquakes? 

None, 

24.  Would  you  be  concerned  if  a major  portion  of  a dam  went  plastic 
under  earthquake  loading? 

No,  but  one  guideline  for  failure  is  if  the  vertical  deformation 
exceeds  0.1  freeboard. 
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Liquefaction 

25.  What  is  the  minimum  earthquake  acceleration  that  can  cause  liquefac- 
tion of  sands? 

Shaking  table  tests  have  shown  that  loose  natural  sands  can  liquefy 
with  .07  g acceleration. 

26.  What  relative  density  is  required  to  withstand  liquefaction? 

Somewhere  around  70  percent  relative  density,  the  material  should  be 
stable. 

27.  How  would  you  assess  relative  density  of  natural  sand  deposits? 

Penetration  resistance  values  are  the  best  available  at  the  present 
time.  However,  there  is  a need  for  a better  technique  to  assess 
relative  density. 

28.  Would  sands  under  a slope  (subject  to  high  shear  stresses)  be  less 
likely  to  liquefy  than  sands  under  level  ground  surface? 

This  answer  may  be  yes,  but  more  work  needs  to  be  done  for 
substantiation, 

29.  Is  there  a decrease  in  susceptibility  of  liquefaction  with  depth  as 
lateral  pressures  increase? 

Yes,  because  the  ratio  of  shear  stress  to  effective  overburden 
pressure  decreases  with  depth. 

30.  What  types  of  laboratory  tests  are  best  suited  to  evaluate  liquefac- 
tion susceptibility  of  sands? 

Shaking  table  tests  with  large  specimens  (2  to  3 ft)  that  are  instru- 
mented for  pore  pressure  measurements. 

31.  What  influence  does  permeability  have  on  progress  of  liquefaction? 

This  has  a large  effect;  with  fine  sand  as  compared  with  gravel  one 
gets  higher  pore  pressures  during  cycling  and,  therefore,  higher 
susceptibility  to  liquefaction. 

32.  What  types  of  field  tests  are  best  suited  to  evaluate  liquefaction 
susceptibility  of  sands? 

No  conment. 
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Laboratory  and  Field  Testing 

34.  How  does  one  interpret  cyclic  load  tests  and  apply  them  in  design? 

No  corment. 

35.  How  does  one  extract  modulus  and  damping  from  cyclic  tests  (if 
possible)? 

No  comment. 

36.  Are  full-scale  field  tests  desirable  and  necessary  for  determining 
dynamic  material  properties  or  are  laboratory  tests  enough? 

These  are  desirable  but  damping  cannot  be  measured  in  the  field. 

37.  For  field  testing,  how  deep  should  one  test  in  a homogeneous  dam? 

(Or  can  one  use  a portable  vibrator  and  save  on  the  shipping  costs 
of  a large  vibrator?) 

Perhaps  one  could  use  a small  vibrator  but  the  effect  of  depth  would 
have  to  be  taken  into  account. 

38.  Is  the  value  of  damping  for  a soil  layer  different  than  that  for  a 
soil  structure? 

Yes,  from  the  viewpoint  of  radiation  damping. 

39.  What  is  the  current  practice  and  application  for  laboratory  tests 
with  respect  to  earthquake  analysis? 

Repeated  load  tests  are  veiy  useful. 

40.  What  are  the  criteria  for  estimating  pore  pressures  in  pervious  shells 
during  seismic  excitation? 

No  comment. 
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Warm  Springs  Dam 


41.  What  earthquake  input  should  WES  use  for  Warm  Springs  Dam? 

An  historical  record,  if  possible,  moved  along  the  appropriate  fault 

to  the  point  nearest  the  site. 

42.  What  are  your  comments  to  the  proposed  analysis  method(s)  for  Warm 

Springs  Dam? 

a.  Determine  susceptibility  of  foundation  materials  to  pore  pressure 
buildup. 

b.  Estimate  amplification  in  structures  and  determine  the  accelera- 
tion levels  that  would  be  present. 

c.  Use  some  technique  to  determine  the  factor  of  safety  along  a 
failure  plane.  If  the  FS  < 1,  use  Newmark's  equations  to  esti- 
mate deformation. 

I 

■ d.  If  a nonlinear  analysis  is  used,  Ambraseys'  or  Newmark's  methods 

are  not  needed. 
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